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We investigate the time-dependent nonlinear optical absorption of a clay dispersion (Laponite) in 
organic dye (Rhodamine B) water solution displaying liquid-arrested state transition. Specifically, 
we determine the characteristic time to of the nonlinear susceptibility build-up due as to the Soret 
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effect. By comparing to with the relaxation time provided by standard dynamic light scattering 
measurements we report on the decoupling of the two collective diffusion times at the two very 
different length scales during the aging of the out-of-equilibrium system. With this demonstration 
experiment we also show the potentiality of nonlinear optics measurements in the study of the late 

o : 

i stage of arrest in soft materials 



■ Light can influence matter microscopic properties providing an optical response that depends on the laser power. 
The nonlinear optical susceptibility of a wide variety of materials (e.g. simple liquids magnetic fluids liquid 
[4 1, polymeric thin films 5] and nano-composite layers jfjj]) has been investigated by using different techniques. 
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Time-resolved Z-scan [7[ has been used to distinguish between different nonlinear effects in crystals and in liquids 
and transient holographic gratings [§] allowed to study the thermal diffusion of nanospheres. Nonlinear optics in soft 



material has been explored in diluted colloidal suspensions 9], in metallic nanoparticles 



3 wi 



with different chemical 



stabilizers 11 1, in presence of nanoparticles aggregation [l^] and recently in out-of-equilibrium systems [lil \l\ 



Complex fluids far from thermodynamic equilibrium arc known to display aging 
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16(| , meaning that dynamical 



quantities such as correlations and responses depend on the time spent since the sample preparation (waiting time 
t w ). The same is expected for the dynamics of nonlinear optical self-action effects, like self-phase modulation and 
nonlinear absorption, an issue that so far remains unexplored. These effects can be used to determine the Soret 
coefficient and the collective diffusion coefficient of the colloidal particles moving on length scales much larger than 
particles size. Therefore the presented nonlinear optics experiments can be very useful to characterize dynamical 
properties of colloidal systems, for instance i) to study the dynamics of jammed systems characterized by long range 
interactions and ii) to analyze the dynamical properties of weakly turbid samples, not detectable by standard light 
scattering measurements. 

In this Letter we report on the time-dependent optical nonlinearity of a complex out-of-equilibrium system during 



its dynamics slowing down. The investigated soft-medium is a dye-doped clay (Laponite) dispersed in water. The 
suspension is prepared by solving the Laponite powder, supplied by Laporte Ltd., in a solution of Rhodamine-B 
(RhB) and deionized water at 0.1 mM concentration. The suspension at 1.6 wt% clay concentration is then stirred 
vigorously until it is cleared and filtered through 0.22 /im pore size Millipore filter. The sample is prepared in air and 
we take the time when the suspension is filtered as the starting aging time (t w = 0). Once dispersed in water Laponite 
platelets start to aggregate and their dynamics slows down leading to a liquid-arrested state transition [lj]]. In the 
mentioned RhB-water-Laponite suspension all dye molecules attach to the clay particles surface We estimate 

about 10 dye molecules per platelet, much less than the maximum allowed adsorption sites (about 500 as provided 
by Laporte Ltd. data-sheet). The RhB molecules absorb visible light that locally heats up the medium and the 
induced thermal gradient influences both water density and dielectric properties (thermal effect) [r^ and colloidal 
particle concentration (thermo-diffusive or Soret effect) 2fJ. Both phenomena lead to a nonlinear optical response 



characterized by distinct time scales: ~ ms and ~ s, respectively. The thermal effect is responsible of local refraction 
changes (thermal lens), while the Soret effect in principle can produce either absorption or refraction variations. To 
study this twofold effect we perform both non- linear refraction and absorption measurements. 

By using the Z-scan approach, for refraction experiments the light source is a CW pumped diode laser operating at 
wavelength A = 532 nm modulated by a mechanical shutter with opening time 700/xs. The beam is focused by a 75 
mm focal-length lens providing a 35 /im beam waist radius w at the sample position. A photo-detector with rise-time 
14 ns and an angular acceptance of 0.73° was used to probe the light power. Being the beam spot on the detector 
larger than the photo-diode surface, only the light transmitted along the beam axis is collected, giving access to the 
nonlinear phase shift 



2l| . For absorption experiments we use the same setup in addition with a 25 mm focal- length 
lens before the detector, thus collecting the whole transmitted beam. The colloidal dispersion is syringed in a glass 
cuvette with transverse dimensions 10 x 35 mm 2 and 1 mm thickness along the beam axis. Each single measurement 
is the response of a different point of the sample not previously illuminated. We investigate the doped clay response 
at four different incident laser powers for each t w , during which we consider unchanged the dynamics and structural 
properties of Laponite at 1.6% w/w, as inferred from Dynamic Light Scattering (DLS) measurements on the same 
sample. 

Refraction measurements are aimed to distinguish between thermal and Soret time-scales. In Fig. [1] the normalized 
transmitted light intensity is reported displaying the two distinct times. The intensity is normalized such that its value 
gives Ah(t) / Aht , where An(t) is the nonlinear phase shift per unit length and unit intensity (radial integral of the 
bell shaped An(r,t), with r the transversal radial coordinate) and Aut corresponds to the saturation of the thermal 
contribution. At short times the temperature gradient induces mostly water density variation, which provides a change 
of the index of refraction A n{t) < (thermal lens) that broadens the input beam. Being the sample positioned before 
the focus of the lens (z=0, as shown in the inset of[T]), the transmitted beam gets narrower augmenting the detected 
light intensity. In Fig. [T] the intensity growth at large times is due to the dyed clay thermal-diffusion, that induces 
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FIG. 1: Refraction. Normalized transmitted light intensity of O.lmM RhB-water+1.6% w/w Laponite dispersion compared 
to the theoretical thermal lens curve (solid line)- see text. Input power 8.3 mW. Inset: sketch of z-scan profile; the cross 
corresponds to the sample position. 



variation of light absorption influencing the refractivity. Z-scans experiments 14] show that the RhB-platclcts diffusion 
drastically alters the nonlinear absorption that is considered with details in the following. In Fig. [TJ at short times, 
the thermal effect is compared to the theoretical curve from 221 ] . calculated with our experimental parameters. From 



the comparison we have t* = 65 ms as the time when the thermal gradient can be regarded as stationary. 

In the following nonlinear optical absorption and its evolution during the gelation of the soft material is considered. 
The transmitted intensity is J(r, t) = Io(r) exp[— a(r, t) L] with Io(r) = Iq exp(— 2r 2 /w 2 ) the input, w the beam waist, 
io the peak value and L the sample thickness. The absorption coefficient a{r,t) = oq + Aa[I(r, t)] can be written as 
the sum of a constant term (the linear absorption coefficient) and a time and laser intensity dependent term Aa[I(r, t)] 
that is related to the gradient of particles mass fraction Ac[I(r, t)] by 



Aa[I(r,t)] 



da 
Ik: 



Ac[I(r, t)], 



(1) 



where c is the average mass fraction. Aa[I(r,t)] is the response of the out-of-equilibrium system to the solicitation 
I(r,t), which is proportional to the local temperature gradient responsible of the platelets diffusion. We consider the 
series expansion of Aa[J(r, t)] in terms of Iq: the first term is the susceptibility linear in theperturbation and it can be 
described by the linear hydrodynamics equations [equations ([2]) and ([3]) below]. The nonlinear absorption coefficient 
Aa[I(r, t)] can be obtained in our experiments by writing Aa(r, t) = —[I(r, t) — I(r, t = t*)]/L I(r, t — t*), where t* is 
the starting time of the colloid Soret effect, I(r, t*) = Io(r) exp[— a(r, t*) L] and we retain the first term in the expansion 
exp[— Aa(r, t) L] . In this way we examine the contribution to the nonlinear optical susceptibility Aa(t) , that is only due 
to platelets motion and hence its dependence on the structural properties of the jelling system. In the following we scale 
down the response function Aa(r,t) to the input power Pq as Aa(r,t) / Pq. In Fig. [2] we show the output signal S(t) = 
— {Aa(r,t)) r /Po at four powers and two different aging times ({Aa(r,t)) r — J Q °°drre ^ 2 Aa(r,t)/ J °° drre w 3 is 
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FIG. 2: Absorption. Scaled absorption coefficient of O.lmM RhB-water+1.6% w/w Laponite dispersion at two different aging 
times. The four inject power values are: 8.3 mW, 11 mW, 14 mW and 17.8 mW. 



the radial average over the incident intensity profile) . In the range 65 ms < t < Is the signals at different input 
powers collapse on one curve defining the linear response regime. For t > 1 s higher orders terms are relevant; these 



23] and their characterization will be reported 



terms give direct access to the dynamical heterogeneity length scales 
elsewhere. 

In the "linear" regime ( 65 ms < t < 1 s) the induced particle mass fraction change Ac[I(r,t)) is given by the 
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solution of the heat-mass flux coupled equations 

^AT(r,t) = D T V 2 AT(r,t) (2) 

D T D D T n c 2 D T a Q 

+ : V Ac r, t) H I(r) 

k k 

^Ac(r,t) = D c V 2 Ac(r,t) (3) 
+ D c S r c(l-c)V 2 AT(r,t) 

where AT(r, t) is the radial temperature variation due to the local light absorption, Dt and k are water thermal 
diffusivity and conductivity respectively, T is the bulk temperature fi c is the the variation of colloid's chemical 
potential in presence of a concentration gradient and Dd is the Dufour coefficient; D c is the particles collective 
diffusion coefficient and St is the Soret coefficient. As shown in Fig. [T]clay concentration gradient builds up once 
the temperature profile reaches a stationary state, this allows to neglect the Dufour effect in the heat equation and 
to take V 2 AT(r, t) in ^ from the stationary solution of @. In decoupling the two effects we consider negligible 
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FIG. 3: Absorption. Scaled absorption coefficient of O.ImM RhB-water+1.6% w/w Laponite dispersion and the fitting curves 
(full lines) by using expression ((Sj . 



the Soret feedback on the temperature profile, being (|ck2 | Io)/oto ~ 10~ D , where the coefficients cxq = 0.86 mm -1 and 
q?2 = —2.8 nm/W are obtained by fitting log [To /lit = Is)] = L [a.® + a 2 To) vs. To. In this approximation the colloidal 
mass fraction variation leads to 



cS T a Po 

Ac(r, t) = 



- Ei\ - 



2r 2 I 



(C 2 



W 2 1 + 2VTD, 

being P = lo^w 2 the beam power, td = w 2 /4D c the characteristic diffusion time over the beam spot size, and 
Ei the exponential integral function. To analyze the observed nonlinear absorption in our experiments by means of 
equation Q we consider the output signal 

S(t) = -{Aa(r,t)) r /P (5) 

= £ (^) 8 ^'-( 1+ 4) 

as obtained by using expressions ((TJ) and ((H). We use equation ((5J to fit the experimental results at different ^ values. 
We show in Fig. [3] the output data with the fitting curves (solid lines) at two aging times and in Fig. |3]-a the fitted 
characteristic thermo-diffusion times td at various t w . We find that this characteristic time exponentially grows 
during aging as the structural relaxation time tj^j (t w ) obtained from DLS measurements, which are contemporarily 
performed on the same sample. DLS data were fitted assuming a correlation function made by square of the sum 
of an exponential function with relaxation time T\ and a stretched exponential function with relaxation time Ti and 
stretching coefficient f3. The mean relaxation time is defined as tm — T2 P^ 1 T (/3 -1 ) with T(x) the usual Euler 
gamma function. Indeed we fit both To(t w ) and the mean relaxation time TM(t w ) (representing the slow dynamics of 
the correlation function) with the exponential growth law r = tq + n exp(fit w ) (full lines trough data) as shown in 
Fig. 0Ja and Fig. HJ-b respectively. Both td and tm describe the collective diffusion of interacting particles on different 
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FIG. 4: Characteristic thermo-diffusion time (a) and mean relaxation time (b) vs. t w and the fitting exponential curves (full 
lines). Inset: the ratio D c /Dm vs. t w 



length scales: DLS probes dynamics on 1/q ~ 50 nm and nonlinear susceptibility onui~ 35 /im. Although we find 
comparable growth rates fj, from the fitting procedure, the ratio D c /Dm [with D c = w 2 /(4td) and Dm = l/(q 2 tm)] 
reported in the inset of Hj-b as a function of t w increases of approximately one order of magnitude in the accessible 
waiting-time window. Moreover at very low t w the two diffusion coefficients are not the same. This behavior can be 
attributed to the long ranged (screened electrostatic) interactions between Laponite particles, characterized by the 
peculiarity of non- homogeneous surface charge distribution, that makes the system strongly interacting on relatively 
large length scales even at very small t u 



251 ] . Another possible explanation might be the deviation to the Fluctuation- 
Dissipation Relation (FDR) in an out-of-equilibrium glassy system [15j. Indeed the diffusion constant measured in 
optical absorption is related to the response of density to an external field, while the one measured in DLS is related 
to the correlation of density fluctuations. The mentioned discrepancy could be also related to a supradiffusive regime 
as shown in 26 1. 

We finally estimate the Soret coefficient as St — 0.02 if -1 from the amplitude of the fitting expression ([5]), by 
considering (^)e — «o/c, and by using water thermal conductivity k = 0.58 W/ m K. Notably enough the Soret 
coefficient, defined as St = — (1/c) d[Ac(t w )]/dT, is constant during the aging process and St > thus indicating 
that the platelets migrate from hot to cold regions being d(Ac(t w ))/dT < 0. 

In conclusion we measured the time-dependent nonlinear optical susceptibility of a soft-material evolving toward 
structural arrest. We characterized the dynamics of the nonlinear absorption in terms of material parameters as the 
characteristic diffusion time td and the Soret coefficient. We have evidence of decoupling of two different diffusion 
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times probed on two distinct length scales signa ling the existence of long ranged interactions between the colloidal 
particles and suggesting the violation of the FDR . 27} . The nonlinear optical response can be hence used to investigate 
dynamical properties of colloidal gels characterized by long range interactions, especially the presented method can 
cover length ranges much larger than particles size, not detectable by scattering measurements; moreover it can be 
really useful in analyzing the dynamics of weakly turbid samples, not easily investigated by standard light scattering 
techniques. 

We thank Giovanni Bolle for helping us in the experimental set-up realization and Francesco Sciortino for useful 
discussions. C.C. acknowledges the ERC Grant No. 201766. 
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